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Abstract 

The general relativistic perturbations of scalar-tensor theories (STT) of gravity are studied in a mani¬ 
festly gauge invariant Hamiltonian formalism. After the derivation of the Hamiltonian equations of motion 
in this framework, the gauge invariant formalism is used to compute the evolution equations of linear per¬ 
turbations around a general relativistic spacetime background in the Jordan frame. These equations are 
then specialized to the case of a flat FRW cosmological background. Furthermore, the equivalence between 
the Jordan frame and the Einstein frame of STT in the manifestly gauge invariant Hamiltonian formalism 
is analyzed, and it is shown that also in this framework they can be related by a conformal transformation. 
Finally, the obtained evolution equations for the linear perturbations in our formalism are compared with 
those in the standard cosmological perturbation theory. It turns out that the perturbation equations in 
the two different formalisms coincide with each other in a suitable limit. 
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1 Introduction 


Perturbation theories of gravity have been very broadly used in both theoretical and experimental explo¬ 
rations during the past few decades. In the theoretical research, linear order perturbation theories have 
been frequently applied in analyzing the stability issues of different kinds of cosmological solutions and 
black holes 0,0, and perturbation theories on FRW cosmological background play an important role in 
explaining the small inhomogeneities of the cosmological microwave background which we have observed 
today Combined with the inflation hypothesis, cosmological perturbation theories enable us to gain 

a possible deep insight into the mysterious early universe 0. 

So far most of the perturbation analysis were done in the context of general relativity (GR). From the 
viewpoints of the solar system experiments, GR is currently the best candidate theory of gravity. Neverthe¬ 
less, for a variety of reasons 0, alternative theories of gravity never cease to exist even from the day GR was 
proposed and attracted more and more attention in the last twenty years ever since the anomalous galaxy 
rotation curves and cosmological acceleration were observed 0, 0 . There are usually two ways to explain 
these two phenomenons. One is to introduce some unknown “dark matter” and “dark energy” components 
which interact with other fields only through gravity The other is to modify the gravity theory. 

While to add a cosmological constant into GR could explain the cosmological acceleration, the extremely 


tiny value of this constant suffers from the “fine tuning problem” 1^. To avoid this problem and also due 


to some other considerations, many people turned to look for modified gravity theories 0. Among all the 
reasonable attempts, scalar-tensor theories (STT) of gravity are the ones that receive most attention. On 
one hand, STT provide t he g reat possibility to solve both the anomalous rotation curve and cosmological 
acceleration problems 0,0. On the other hand, STT can include a lot of modified gravity models as its 
special sectors, such as f{R) theory, Brans-Dicke theory, induced gravity, etc 0,[i3 . What is worth men¬ 
tioning is that, when studying the cosmological perturbations during the slow-roll inflation period of STT, 
people noticed that the physics there behaves much different from that of GR 0, 0]. Some predictions 
made by f{R) = R + aR^ gravity and certain STT theories can predict the scalar spectral index and the 
tensor-to-scalar ratio completely consistent with the Planck data [20|. A question that cannot be answered 
in the context of STT is how matter couples to gravity, but is usually treated as an additional freedom of 
the model under consideration. The two most prominent choices for STT are the Einstein and the Jordan 
frame, which are related by a conformal transformation, and which will be also chosen in this paper. A 
systematic approach how alternative actions for gravitational theories can be derived once the information 
how matter couples to gravity is known, can for instance be found in 2ll |. 


GR can be understood as a gauge theory with the gauge group being the diffeomorphisms Diff{M) of 
the spacetime manifold. Likewise to other gauge theories physically relevant quantities, called observables, 
are those that are gauge invariant, and in the context of GR this corresponds to diffeomorphism invariant 
quantities. In the canonical formalism, that is used in this paper, this carries over to the requirement that 
observables, also denoted as Dirac observables in this context, are tensors on phase space that do (weakly) 
Poisson commute with all constraints of the system, and particularly for GR in the ADM formalism this 
means, with the Hamiltonian and spatial diffeomorphism constraints. The elementary variables in the ADM 
formalism are the ADM 3-metric, given by the pull back of the spacetime metric onto the spatial hyper¬ 
surfaces, and its conjugate momenta. The canonical Hamiltonian density consists of the sum of the spatial 
diffeomorphism and Hamiltonian constraints, weighted by the so called shift vector and lapse function re¬ 
spectively, which parametrize possible foliations into spatial hypersurfaces of the four dimensional spacetime 
manifold. These configuration and momentum variables are no observables, but their Poisson brackets with 
the canonical Hamiltonian, and hence with the constraints, together with the constraints themselves yield 
the analogue of Einstein’s equations in the canonical framework. As a consequence, when we discuss ob¬ 
servables in the context of GR, the natural question arises, how the evolution of those observables can be 
described. Gertainly, it cannot be generated by the canonical Hamiltonian, since this would allow only trivial 
evolution contradicting what we observe, because we do observe the evolution of physical observables in our 
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everyday life and in experiments. Thus, what we are looking for is a so called physical Hamiltonian having 
the property that it generates the evolution of the observables with respect to physical time in the canonical 
setup of GR. The initially missing physical time and along with this the missing physical Hamiltonian is a 
common feature of any diffeomorphism-invariant theory and is often called the ’’problem of time” for short. 
Many attempts have been made to clarify this conceptual issue in the history, see for some examples. 
One possible way to circumvent the problem of time was introduced by Rovelli and mathematically fur¬ 
ther developed by Dittrich 2^ in the context of so called relational observables. The main idea of relational 
observables is to introduce so called reference fields with respect to which the observables and the evolution 
of the remaining degrees of freedom are formulated. These reference fields could for instance be additional 
matter fields but also purely gravitational degrees of freedom. In a seminal paper by Brown and Kuchaf 
they considered to use pressureless dust particles as such reference fields, which were introduced as additional 
matter fields. The dust particles can be understood as a free falling observer that is dynamically coupled to 
the system and with respect to which observables for GR are constructed. The dust fields therefore serve 
as a physical reference system and provide an interpretation of physical spatial and time coordinates. In 


25l | observables were constructed explicitly only with respect to the spatial diffeomorphism constraint and 


for the Hamiltonian constraint a corresponding Schrodinger-like equation was discussed within the quantum 
theory. A combination of the methods of Brown and Kuchaf and the relational observable framework was 
used in 2^ to derive a manifestly gauge invariant Hamiltonian formulation of GR. The corresponding phys¬ 
ical Hamiltonian associated with the dust observer then generates evolution with respect to this particular 
chosen physical time. As an application in 2^ the manifestly gauge invariant Hamiltonian formalism was 
used in order to develop a formalism for general relativistic perturbation theory. The obtained results were 
then in 2^ applied to the context of cosmological perturbation theory for which a comparison with the 
results of the standard formalism, in the following denoted as standard cosmological perturbation theory 
(SGPT), was analyzed for linear order perturbation theory. It turns out that they are consistent with each 
other only up to small corrections. The latter occur due to the reason that in the manifestly gauge invari¬ 
ant formalism the observer is dynamically coupled to the system and thus has an influence on the system, 
whereas this is not the case in SGPT. 

The main difference between the manifestly gauge invariant and the standard formalism lies in the way how 
gauge invariant quantities and thus observables are constructed. In the manifestly gauge invariant formal¬ 
ism first one constructs observables and the associated physical Hamiltonian, which itself is an observable. 
Then one derives the equations of motion for all observables and afterwards considers the perturbation of 
these evolution equations. By construction any quantity and their corresponding perturbations, which are 
involved in these evolution equations are themselves observables and hence gauge invariant. Following SGPT 
the strategy is different since perturbations of gauge variant quantities, the metric and and their equations 
of motion, are considered. Gange invariant objects, and hence observables, are constructed afterwards and 
are required to be only gauge invariant up to the corresponding order of perturbation theory one is inter¬ 
ested in. The manifestly gange invariant formalism might be of advantage when higher order perturbation 
theory is considered because in SGPT one has to start from scratch again when gange invariant quantities 
are constructed [^, while using the manifestly gauge invariant framework any object and thus also higher 
order perturbations are already gauge invariant. 


In this paper we extend this manifestly gange invariant formalism to the case of STT. This paper is organized 
as follows. In section two, after introducing the Brown-Kuchaf formalism, relational observables and the 
notion of a physical Hamiltonian, we construct the physical Hamiltonian in Jordan frame of STT and derive 
the second order evolntion equations of the canonical variables. In section three the latter will be used to 
derive the evolution equations for linear perturbations in the context of a general relativistic spacetime. In 
section four we apply these equations to flat FRW background and get the cosmological pertnrbation equa¬ 
tions. In the first part of section five we formulate the gauge invariant Hamiltonian formalism in Einstein 
frame of STT and prove it is conformally equivalent to that of Jordan frame. In the second part we extend 
our results and stndy the cosmological perturbations in a different reference system. Our results are then 
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compared with the standard cosmological perturbation theory of STT. In the last section we summarize the 
above results and draw some conclusions. In the appendix the Lagrangian and Hamiltonian analysis of the 
action of STT in Einstein frame is presented. 


2 Manifestly Gauge Invariant Hamiltonian Formalism of STT 


For the convenience of readers, in the first part of this section, we will first give a brief review of the idea pro¬ 
posed by Brown and Kuchaf. Then we introduce how to combine the Brown-Kuchaf’s dust formalism with 
the idea of relational observables to build the physical reference system and the gauge invariant Hamiltonian 
which generates the evolution of every observable in this physical reference system. We refer to 2g, [29| for 
the details. In the second part we will apply these ideas to STT of gravity. First we derive the physical 
Hamiltonian in the Jordan frame of STT and then use it to derive the evolution equations of the 3-metric 
and gravitational scalar field. 


2.1 The Brown - Kuchaf Lagrangian and Relational Observables 

In this subsection we will briefly review the idea of Brown and Kuchar to use dust as reference system for 
general relativity and discuss how this idea can be embedded in the framework of relational observables. 
The generalization from general relativity to STT is not difficult and will be discussed in section 12.31 


2.1.1 Brown — Kuchaf Lagrangian and Deparametrization of General Relativity 


In 2^ Brown and Kuchaf followed the idea that matter fields can be chosen as a physical reference systems. 


For the reason that they were mainly interested in GR they considered the Einstein-Hilbert action and the 
following additional matter action 


Sdust = -\j + 1) 

M 


( 1 ) 


which, as we will see below, can be interpreted as an action for pressureless dust. The action S^ust is not taken 
as a functional of the one form U but the latter is expressed in terms of the scalar fields T, , Wj defined 
through U = —dT WjdS^ where we use the notation that Latin letters run from 1 to 3 and Greek letters 
from 0 to 4. Hence, the action above is a functional of the fields p, T, , Wj and g^i, and hence in addition 
to GR we have introduced eight more degrees of freedom. As we will see below the system GRj-dust has 
second class constraints and once these are reduced the additional number of degrees of freedom is reduced 
to four. As discussed in detail in [25(] the Euler-Lagrange equations associated with this action show that 
the vector field = g^'^Uy satisfies the differential equation of a geodesic in affine parametrization. The 
fields Wj and are constant along the geodesics and the field T defines proper time along each geodesic. 
Furthermore, the energy-momentum tensor associated with S^ust has the form of the energy-momentum 
tensor of a pressureless perfect fluid. In order to use the dust matter as a dynamically coupled observer 
in the canonical framework we have to discuss the Hamiltonian formulation of which is discussed in 
detail in ( 2 ^ . We assume that M is globally hyperbolic and thus we can perform a 3-|-l-split of M ~ M x y 
into time and space, with y being a spatial manifold of arbitrary topology. For this purpose we introduce 
a family of embeddings Xt : x ^ ^ x ^ Xt{x) := X{t,x), also called a foliation of M, where xt '■= ^t{x) 
are called the leaves of the foliation, where we denoted the coordinates on y by with a = 1,2, 3. Given 
the family of embeddings Xt we can construct a family of tangent vectors and a co-normal for each 

leave. Using the metric we can also work with the future orientated normal and can use it to decompose 
the variation of the embeddings with respect to the parameter t into a tangential and normal part to the 
leaves xt given by <9*Af = nn^ -|- n“Af^ where this decomposition is parametrized by the so called lapse 
function n and shift vector The spatial three metric qah intrinsic to x can be constructed by pulling 
back g^j_iy using the tangent vectors yielding qab = g^i^X^aX''^. In order to derive the Hamiltonian formulation 
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via the Legendre transformation we introduce conjugate momenta p°‘^,p,Pa,P,Pj,I,P associated with the 
configuration variables Qab, n, n“, T, , p, Wj. Note that in case we would like to couple additional standard 
model matter to gravity we would need to introduce additional phase space variables for this matter, which 
we will not display explicitly here, but we will mention below how such additional matter degrees of freedom 
enter the model. The theory possesses the following primary constraints 


z ■= p = 0^ Za := pa = 0, Z := I = 0, := P = 0, 


Zr-=Pj-PWj=0. 


(2) 


Following the Dirac procedure for constrained systems the stability analysis of the primary constraints with 
respect to the primary Hamiltonian yields the following secondary constraints 




^ _j_ ^dust ^dust _ ^ 


C = Ca + cf = P[Ta - 


p2 


pVdet{q) 


+ pv^det(g)(l + q'^’^UaUb) 


n 


c = 


p2 


pVdet(g) 


+ v'det(g)(l + q'^’^UaUb) 


(3) 


here c denotes the gravitational contribution as well as the contribution from other possible standard model 
matter to the usual Hamiltonian constraint and Ca denotes the contributions from gravity as well as other 
possible standard model matter to the usual spatial diffeomorphism constraint. Furthermore we have Ua = 
—T^a + WjS^a- A further application of Dirac’s constraint algorithm shows that no tertiary constraints occur 
and thus the set of constraints has completely been determined and is given by {c*°*, c, z, Za, Z^ Z^^ Zj}. 

The next step is to classify them into first and second class constraints. The constraints {c, Z, Z^, Zj'\ form 
second class constraints. We will solve them explicitly below using the associated Dirac bracket. Furthermore 
we will also solve the constraints z,Za by considering the lapse function n and the shift vector n“ as Lagrange 
multipliers. Now we extend and to and by adding terms proportional to the constraints 
z, Za, Z, Zp Zj in the case of and terms proportional to Z,Zj,Z^ j^the case of c*°L These terms are 


exactly chosen such that and are first class constraints, see |26|] for more details. Since Za,z are 
hrst class constraints as well we end up with the following set of first class constraints {C^°^z, Za}- 
solving the second class constraints strongly leads to 


Now 


Wj:=^, p^:=^(^q-^UaUb + l 


I ■= 0, P := 0. 


(4) 


For the reason that Z, Z^ and Zj do only depend on the dust variables the associated Dirac bracket reduces to 
the Poisson bracket when applied to the geometrical degrees of freedom and other possible standard 

model degrees of freedom. When we introduce the Dirac bracket and also solve the constraints z, Za by 
considering the lapse function and shift vector as Lagrange multipliers, as usually done the ADM framework, 
we work on a reduced phase space where and Inserting the explicit solutions of the 

second class constraints from dl]) we end up with the following first class constraints 




Aust 


= P‘^ + q' 


ah r>dust rAust 


J.ot 


= Ca + d. 


dust 


Aust 


= PTa + PjSla- 


(5) 

( 6 ) 


Note that in principle we have two possible choices for the sign of p here but the one chosen by us p < 0 
yielding also to P < 0 is the part of the phase space that involves also flat space solutions since for them 
c > 0 is necessary, for a more detailed discussion about this aspect see (26|. 


The Hamiltonian and spatial diffeomorphism constraints satisfy a complicated Poisson algebra also called 
the hypersurface deformation algebra for the reason that it can also be derived from purely geometrical con¬ 


siderations in the context of deformation of hypersurfaces 30|] • One of the motivations in 25|] to introduce 
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the dust as a reference systems was that using the dust one can write down an equivalent set of first class 
constraints that has the property that the corresponding constraint algebra becomes Abelian and the final 
Hamiltonian constraint can be written in deparametrized form, as we will discuss below. The important 
observation by Brown and Kuchaf, also denoted as the Brown-Kuchaf-mechanism in 2y, [271], was that on 


the constraint surface of the spatial diffeomorphism constraint = 0, we have and thus we 

can rewrite the dust contribution as 


+ (7) 

As a consequence the only dependence in on the dust variables is via the dust momentum P. Thus on 
the constraint surface of the Hamiltonian constraint we can solve for P and obtain an equivalent form 
of the Hamiltonian constraint given by 

= P + h, h:= i/c2 - q°-’^CaCb. (8) 


The fact that h no longer depends on the dust variables is what is called deparametrization and a consequence 
of this is that we will end up with a time independent physical Hamiltonian as discussed in the next 
subsection. In order to obtain an Abelian constraint algebra we also solve the spatial diffeomorphism 
constraint for the momenta Pj. For this purpose we have to assume that the matrix S^a is everywhere non¬ 
degenerate, an assumption similar to the classical restriction det(g) > 0, meaning that the inverse matrix 
exist, which we denote by Then on the constraint surface we have Pj + Sj{ca + PT^a) = Pj + Sj{ca — hT^a) 
thus we can write down the following equivalent form of the spatial diffeomorphism constraint 


-tot 


— P-j h 




hj := - hTa). 


(9) 


We realize in contrast to the Hamiltonian constraint the spatial diffeomorphism constraint does not 
deparametrize. However this is no problem at all because the construction of observables is not restricted to 
the deparametrized case and can therefore be equally well applied to the spatial diffeomorphism constraint. 
It might only be technically a little bit more involved. It is only at the level of the physical Hamiltonian, 
for which however only the form of turns out to be important, where deparametrization yields to 
simplifications in the sense that the final physical Hamiltonian will be time independent. 

Now considering the constraints and one can indeed show that they satisfy an Abelian constraint 


algebra 25|. This follows also immediately from the following abstract argument 3l||: The equivalent 


constraints and are still hrst class. Therefore their Poisson brackets are again linear combinations 
of constraints. However since all constraints of the system are linear in the momenta P, Pj their Poisson 
brackets are independent of P,Pj- Consequently, we can evaluate the linear combination of constraints that 


appear in the Poisson bracket computation in particular at P = —h and Pj 


— hj. 


From the Abelian 


constraint algebra and the explicit form of h and hj in ([8]) and ([9]) respectively we can conclude that h{x) 
are mutually Poisson commuting while h{x) does not Poisson commute with hj{x) nor do the hj{x) mutually 
commute. In the next subsection we will introduce the relational framework and use the latter to construct 
observables in the context of GR. 


2.2 Relational Framework for Constructing Observables 

The main idea of the relational framework is to introduce so called reference fields, often also denoted as 
clocks, that will then be used to construct observables with respect to the constraints of the system under 
consideration. Let us assume we have a system with a set of constraints {C/} labeled by an index /, which is 
up to now arbitrary. The aim is to introduce for each constraint Cj a corresponding reference held such 
that the constraint and the reference held build, at least weakly, a conjugate pair, that is {C' 7 ,r‘^} « 
where ~ means equality up to terms that vanish on the constraint surface. Now since for a given set of 
constraints hnding those reference helds might not be a simple task, one uses the freedom that one can 
always modify the set of constraints as long as the modihed set dehnes the same constraint surface. Suppose 
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we choose a set of reference field {T^}, one for each Cj with the property that {Ci,T'^} =: Mj with M 
being an invertible matrix, then we can define the equivalent set of constraints {Cj} defined through 

(10) 

J 

One can easily show that for {C'j} we have {C'j,T'^} ss 5j. Given these new set of constraints {C'j} we can 
use the reference fields {T^} to construct observables for a general phase space functions. This will be a 
particular combination of the original phase space function under considerations and the reference fields. To 
discuss this construction more in detail we consider the Hamiltonian vector field associated with Cj, which is 
denoted by Xj. As can be shown and will be crucial in the following constructions the Xj mutually weakly 
commute. Let us introduce a set of up to now arbitrary real numbers {/3/}, again one for each constraint 
Cj, and consider the following sum of Hamiltonian vector fields 

Xp-.= Y,P'Xi. (11) 

I 

Now we consider a function / on phase space and define a map / «/?(/) on the set of smooth functions 

on phase space given by 

CXD . 

a0{f) := exp(X^) . / = ^ • / , (12) 

1 

ri.=l 

here X^ ■ f = {Cp, /}(n) where .}(„) denotes the iterative Poisson bracket defined through {Cp, /}(o) = / 
and {Cp,f}(^n) = {Cp, {Cp, f}■ cup is a Poisson automorphism on the algebra of functions on phase 
space associated with the Hamiltonian vector field Xp of Cp = P^Cj. We will use the map ap as well as 
the set of reference fields to construct an observable associated with a phase space function /. A weak 
Dirac observable has to weakly Poisson commute with all constraints {G/}. Now the idea of the relational 
observables is that although the phase space function / as well as the reference fields have non-vanishing 
Poisson brackets with the constraints a particular combination of the two involving the map ap has vanishing 
Poisson brackets with all constraints. We want to construct a map that returns the value of / at those values 
where the reference fields take the values r^. In order to do so let us choose another set of real numbers 
{r^}. We are interested in those values of the gauge parameters for which ap{T^) = r/. If we apply ap 
onto the reference fields we obtain ap{T^) k,T^ + (3^, which can easily be solved for f3j yielding fjj = . 

We will denote this equation for short as (5 = t — T suppressing the indices. Using this we can construct 
the following map for the phase space function / 


0 /( t ) := [ap{f)]p^^_r^ 


(13) 


The notation with the square brackets means that only after one has computed the action of Xp with j3 
treated as a constant on phase space then one sets f3 = t — T which becomes then phase space dependent. 
As has been proven in 2J, [3^ Of{T) is indeed a weak Dirac observable, that is for all I we have 




(14) 


We realize that we can also understand the map Of as a map that returns the value of / in the gauge 
P = t — T. As also shown in 3l|, [3^ the multi parameter family of maps : f Of^r) is a homomorphism 
from the commutative algebra of functions on phase space to the commutative algebra of weak Dirac 
observables, both with pointwise multiplication. 


0/(r) -h Ogir) = 0/+g(r), Of{T)Og{T) « 0/g(r). 


(15) 


This will be a particularly useful property when the explicit construction of the observables is considered 
for the following reason: Let us denote the coordinates on phase space by {q^,PA) where the index A is 








chosen such that all relevant phase space degrees of freedom are considered. Now for a phase space function 
/ = have we have 

Ofir) = f{Oq^,OpA){T) (16) 


which has the important consequence that it is sufficient to construct observables for the elementary phase 
space variables, something we will use below. Moreover, multi parameter family of maps O'^ : f -A 0/(t) 
is a Poisson homomorphism with respect to the Dirac bracket {., .}* associated with the system of second 


class constraints Ci,T^ 3^ , this means 


{0/(t), Og{T)} « {0/(t), Og{T)Y K (r) (17) 

where the Dirac bracket is defined as 


{f,gr = {f,g} - {f,Ci}{M-Y^j{TYg} + {g,Cj}{M-Y^j{TYf}. (18) 

In the following we want to discuss the special case of constraints that are in deparametrized form and 
understand how this simplifies the construction of the observables Oj(r). In the case of deparametrization we 
can always find canonical coordinates that consists of two sets {T^,Pj) and (g“. Pa) such that all constraints 
Cj of the system can be written in the following form 


Cl = Pi + hi{q°‘,Pa), (19) 

and thus do not depend on the configuration variables . In practice this is a very special case and most 
constrained systems, if at all, can only be written in partially deparametrized form, in which only part 
of the constraints deparametrize. However, for the following discussion let us assume that we consider a 
fully deparametrized system. Now following the steps of the construction of observables from the discussion 
above first we observe that 

{Ci,TY = 6j. ( 20 ) 

Using the notation above this means the equivalent constraints C) are identical to Ci and thus the task 
of inverting a in general complicated matrix M'/ is no longer necessary. Furthermore as already discussed 
above if all constraints are linearly in the momenta Pi then the associated constraint algebra is Abelian. For 
the reason that here also non of the hi depends on the reference fields T/ we immediately get {hi, hj{ = 0 
from this we can follow {/i/,C'j} = 0 showing that each hi is already a Dirac observable. Moreover from 
the Abelian constraint algebra it follows that also the associated Hamiltonian vector fields commute and in 
this case here not only on the constraint surface but on the entire phase space. As a consequence all weak 
equalities that we used above can be replaced by strong equalities here. 


First let us discuss the construction of the observables for the elementary variables {q°',Pa)- Since and 
Pa both commute with all momenta Pj we can consider the Hamiltonian vector field associated with the 
hi’s instead of defining via C'j. Moreover for the reason that also q^" and pa commute with all reference 
fields we can already, when applying Xp to /, replace /3 by the corresponding gauge Y — yielding the 
following form for the observables for a function / that depends only on {q°‘,Pa) 

OO ^ 

(21) 

n=0 

where Xj- is the Hamiltonian vector field of the function 

Hr = (r^ - T^Hi (22) 

where Hi := Ohjir) denotes the observables associated with hi. Because hi = hi{q^,pa) is a function of 
q°' and pa only, once the observables for the elementary variables Oq^Y) =: Q“(r) and Op^Y) =: PaY) 


9 




are constructed we obtain Hj as Hj = Ohj{T) = hj{Q°’, Pa){T) using the homomorphism property of the 
observable map. In the particular case of deparametrization we have Hj = hj because hj is already a Dirac 
observable as discussed above. Now if we restrict to functions that do only depend on and Pa the Dirac 
bracket reduces to the Poisson bracket because those / commute with all reference fields . In particular 
for the algebra of the observables Q'^{t) and Pai^) we obtain 

{Pa{T),Q\T)} = {OpJr),0,.(r)} = = 0,.(r) = 5^ , (23) 

showing that the reduced phase space has a very simple symplectic structure in terms of the coordinates 
Q“,Pa, an important property if the quantization of such systems is considered. Having finished the discus¬ 
sion about the non-reference held degrees of freedom let us discuss now the case of the remaining reference 
held degrees of freedom. The observable associated to the reference helds is given by 

Oti{t) = (24) 

and therefore is just a constant function on phase space. Since all momenta Pj Poisson commute with all 
constraints they are already Dirac observables. In addition they can also be expressed as function of the 
observables Q“(r) and Pai^), because on the constraint surface we have 

Pi = Op,{t) = -Ohjir) = -hj{Q^{T), Pair)) = -Hj . (25) 

Hence, what we hnally be interested in is the reduced phase space with elementary variables Q“(t) and Pai^)- 

Let us again consider an observable Of{T) associated with a function that depends only on q°’ and pa- 
How can we formulate the evolution of such observables? Certainly this cannot be generated by the con¬ 
straints since by construction 0 /(t) Poisson commutes with all constraints. However, 0 /(t) gives us the 
value of / when the reference helds take the values r^. As it will be the case for GR and also for STT 
one of the chosen reference helds will be associated with physical time and let us without loss of generality 
denote this reference held by and the values that it takes by Then time evolution for Of{T) can be 
described by the derivative of 0 /(t) with respect to since this encodes how 0/(r) changes with time r®. 
However, considering the form of 0 /( t ) in (I2ip we can explicitly compute this derivative and as shown in 
[ 3 ^ one obtains 

(26) 

OTq 

where PIq is the observable associated with ho that occurs in the constraint Cq := Pq + ho associated with 
the reference held that we interpret as a reference held for time. In the following we will call Hq the 
physical Hamiltonian because in contrast to the constraint Cq, that is generating gauge transformations, 
PIq does not vanish on the constraint surface and can therefore be understood as a true Hamiltonian, which 
generates evolution with respect to physical time r°. Note that because ho does not depend on (and also 
not on any other reference held) the hnal physical Hamiltonian PIq is time independent. 

2.2.1 Observables for GR Using the Brown-Kuchaf Dust 

In this subsection we will discuss how the Brown-Kuchaf dust can be used to construct relational observables. 
In the case of GR and STT we have four times inhnitely many constraints because we have one Hamiltonian 
and three spatial diffeomorphism constraints per spacetime point. Following the discussion of the last 
subsection we therefore need to choose 4 times inhnitely many making four scalar helds a natural choice 
for reference helds. These will become exactly the four additional degrees of freedom (T, S^) which we added 
to the system by considering the Brown-Kuchaf Lagrangian. The remaining degrees of freedom {qabiP^"^) 
and possible other standard model degrees of freedom will be referred to as non-dust degrees of freedom. 
In the following in order to keep the discussion more simple we will only consider the system of gravity and 
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dust. In order to define the Hamiltonian vector field Xp in this case, we introduce arbitrary functions /3^, /3-^ 
on X aiid define using the constraints in ([8]) and ([9]) 

cf := J d^xpi^{x)c^°\x) (27) 

X 

where /3^ = {(3^,(3^) and we have defined = c*°*. We denote the Hamiltonian vector field of by Xp 
and using it we can define the map given by 


- 

«/?(/) := exp{Xp) ■ f = Y^ —X^ • / . (28) 

n=o 

We use the notation r^(x) = (T(x), S-^(x)) then applying onto yields ap{T^^{x)) = T^(x)+/3^(x). Let 
us denote the values that the reference fields T^{x) can take by r^(x) with t°(x) = r(x) and r-^(x) = a^{x) 
where up to know are arbitrary functions on y. Now solving ap{T^{x)) = r^(x) for f3^{x) leads to 
/3^(x) = r^(x) — T^(x), which we again write as j3 = t — T. The observable associated to / reads 


Ofir) = . (29) 

Looking at the explicit form of in ([8]) we realize that c*°*(x) commutes with S^{y) and because of this 
we can construct the observables in two steps. First we reduce with respect to the spatial diffeomorphism 
constraints and afterwards we construct the complete observables that also Poisson commute with the 
Hamiltonian constraint. Hence, we can rewrite 0 /(t) as 


Of{T,a^) 


a^o ((/)] 


I30=t-T 


(30) 


Let us first discuss the inner part, that is how spatially diffeomorphism invariant objects are constructed. The 
reference fields will be used for this and therefore here the remaining degrees of freedom are {qab,P°‘^, T, P) 
for which observables need to be constructed. As discussed in detail in [ 2 ^ for the choice of a constant 
function cr^ (x) = the observable with respect to associated with any scalar function / build from the 
variables {qab,P°‘’^,T, P) can be expressed as 


/V) = 


fix) 




(31) 


where we denote the partially reduced function as /. The interpretation of the formula above is the following: 
Whatever the value x is at which the function / is evaluated /(cr-^) returns / evaluated at the point x^ at 
which S^{x) = aP We call the range of the dust space and denote it by 5. Since by our assumption 
is everywhere invertible it defines a diffeomorphism : x —)■ 5 and hence the value of x^ is unique. Our 
strategy is therefore to use to construct scalars f on x for {qab,p°'^,T, P) and then apply the formula in 
m yielding the partially reduced quantities {qab,P°'^,T, P), explicitly we get 


T = T 




qjk = qabS^St 


(32) 


where J := det (§§) was used to obtain the correct density weight. Note while these are scalars on x there 
are tensors on the dust space S with the same density weight that they have on X; see also 2^ for more 
details. We realize that the evaluation of the functions in (|32p at x^ is nothing else than the pull back of 
the corresponding fields to S under the inverse of the diffeomorphism : x ^ S. 

Our remaining task is to compute the complete observables in (I30p that also Poisson commutes with the 
Hamiltonian constraint. For this purpose the reference field T will be used and hence we need to construct 
observables for qjk and that we will denote by Qjk and P^^ respectively. The constraint cf' 


-tot 


in 


IS m 
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deparametrized form and thus wen can apply the simplified construction discussed in the last subsection. 
Let us look at the smeared version of the constraint given by 

c“ := j d^x{T - T){x)c^°\x) (33) 

X 

where we used = (r — T){x). Now in order to construct the observables Qj^ and we need to 

ensure that is already an observable with respect to the spatial diffeomophism constraint . For the 
choice of constant r the constraint is an integral over a density of weight one and we can equivalently 
express it as an integral over the dust space S given by 

cf = J dV(r - f)ia)iP + h)ia) (34) 

5 


with 

h{a) = ^ c 2 (( t ) - qi^{a)cj{a)ck{(T) ( 35 ) 

where c{a) and Cj{a) are the observables of c and Cj in ([5]) and @ respectively with respect to the spatially 
diffeomorphism constraint. Now since deparametrizes we do not need to consider the Hamiltonian vector 
field of but we can work with the Hamiltonian vector field X-r of 

Hr ’■= J (Pa{T — T)h{a) . ( 36 ) 

5 

The observables for a function / that depends only on qjk and (and possible other standard model 
matter degrees of freedom) is then given by 


oo ^ oo ^ 




(37) 


n=0 


ri.=0 


where {.,.}(„) again denotes the iterated Poisson bracket. Considering the discussion about the physical 
Hamiltonian above we observe that in case of the dust as reference fields the physical time is given by r and 
the corresponding evolution of the observables is given by 


dOf{T,a) 

dr 


{H,0/(r, 0 -)} 


with the physical Hamiltonian 



Here H[a) is the (complete) observable associated to h{a)^ that is 


■= O-hiT.a) = h{a) 


(38) 


(39) 


(40) 


and is therefore independent of r and hence physical time as expected in the deparametrized case. In the 
following we will continue to use the notation Qjk and for the observables of qab and Furthermore we 
will denote the observables associated with c and Cj also by capital letters C'(r, u) and Cj{T,a) respectively. 
Using this notation we can also rewrite the physical Hamiltonian density as 

H{a) = ^C2(r,cr) - Q^pT,a)Ci{T,a)Cj{T,a). (41) 
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2.3 Evolution Equations in Jordan Frame 

As is well known, in the classical formulation of STT, there are different choices for basic variables, corre¬ 
sponding to different frames. These frames are related to each other by conformal transformations. There 
are many debates on which frame should be regarded as the physical one (the word “physical” has nothing 
to do with the gauge invariance here) and the question carries over to the discussion what kind of coupling 
should be chosen. Two prominent choices are either the Jordan frame, where a scalar field is non-minimally 
coupled to the metric, or the Einstein frame, where a scalar field is minimally coupled as in GR. If one 
considers the Einstein frame as the relevant one, the corresponding STT would make no difference from GR 
with a scalar field, which was already presented in [^. Hence, to analyze the non-minimally coupled case, 
in this paper, we consider the Jordan frame and this choice has the consequence that also the evolution 
equations of the linear perturbations will be formulated using the Jordan frame. 

We consider the following action 


SsTT — St Jordan T Sjustj 
where the STT action in Jordan frame reads 


S Jordan = / d^X,/\d^\ - K ,cj) - R(<()) 


(42) 


Here for simplicity we set IGvrG = 1, F{<p) and K{(j)) are positive coupling functions, and V{4>) is the 
potential of the gravitational scalar field (j). As mentioned in section 2.1, the dust action reads 


Sdust = j det( 5 )|/o [g^'^U^Uu + 1 ] 


(43) 


with Ujj = —T/j + WjS^n- The Hamiltonian analysis of STT in the generalized Brans-Dicke form is given in 


Refs. 33| in Jordan frame. Using the results there and following the procedures in the last subsection, the 


physical Hamiltonian density reads 


H{a) = a) - {T,a)Ci{T, cr)Cj{T, a), 


with 


Cj{T,a) = 
C(r,cr) = 


-2QjiDkP^^ (r,(T), 

1 


(44) 


(45) 


(QimQjn - ^QijQmn]P^^P^ 






VdetQV F(4>) ' 2T(4>)[3(F'(4>))2+ 2F($)i^($)] 

+Vdetg( -F(4>)R(3) +i;C($)g*^(A^)(^i$) +2 Q*^AT»j-F(4.) + U($)) 


T, cr 


(46) 


and the elementary Poisson bracket read; 

{P^^{a),Qki{d)} = {n(u), <I>(d)} = J(u,d), (47) 

where a prime ' in Eq. (1461) means the derivative with respect to 4>, i.e. T'(<h) = , F"{^) = Here 

we denote all the objects after gauge completion with capital letters. All quantities are now evaluated on 
the dust manifold. Note that Eq. ([JG)) is valid in the case of 3(T'(<I>))^-|-2T(<h)Ar(<h) / 0, which corresponds 
to the case a; / — | in Refs. (33|, on which we focus in the paper. 
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Given a function / on phase space, that does not depend on the dust variables, we can construct its 
associated observable O/ by applying the procedure described in section 12.11 We get the Hamiltonian 
equations of motion for Of, being first order differential equations with respect to the physical time r, by 
simply calculating its Poisson bracket with the physical Hamiltonian denoted by H and given by H := 
fdaH{a): 

= ^ = (48) 

In our case we are interested in the equations of motion for =: Qjk, Opab =: =: <I>, Ojr =: H and 

using the physical Hamiltonian of the dust model we end up with the following Hamiltonian equations: 


6 

H 


Qjk 


pjk 


N .r ^ 

Vd^3(F'($))2 + 2F($)i^($) ^ N ’ 


(49) 


N 


F'(4>)QijPd' F{^)Uy 


VdetQ L2(F($)) [3(F'($))2 + 2F(4>)K(4>)] 
F'i<^)QijP^^ - F(4>)n 


+ - 


- F(4>)(3(F'(4>))2 + 2 F{^)K{^)) 


F"(4>)Q*jP*^' - F'(4>)n) 


pij pr. 


F(4>) [3(F'(4>))2 + 2F(4>)it'(4>)] 

F'(4>) ^ ^r^ r) 

~ (F(¥))2 ~ 

-y/detQ [-A^F'(4')i?(3) + + 2F\^)Q^WiDjN + NV’{^) 

+2dj 


N^ydet 


N 


VditQF($) 


y^GjkmnP^^ + 


(F'm^QmnP^^ - F'($)F(4>)n 


3(F'(4>))2 + 2F(4')iir(4') 




N 


-v/det Q 


Qmn( 2 P^^P'^^ - 


m) 


(F(ch))2Q^,p^’^ - F(.h)p(ch)n ^^.fe 
P(4>) [3(P'(4>))2 + 2P(d>)P(4>)] 


(50) 

(51) 


+ ^/de^Q 


N 


- Q^^DWiF{^) + 


+ (^D^DniNF{^)) - NF{^)Rmn) - 2{D^^N)D^'^F{^) 


+Q^^{D^N)DiF{^) 

+ ^qFc{t, a) - ip(r, + (£^P)'^ (52) 

where N := -jj, Nj .— -p, 4* ;;; .— Gjkmn ■— 2 (^ 1 ™^^^ 4“ QjnQmk QjkQmn) nnd its inverse 

:= — 2gPQ”*”') satisfying GjkmnG'^'^'^^ = indices are lowered 

and raised by the three-metric and its inverse. Note that both G{t, a) and P(t, a) are non-vanishing because 
they are no longer constraints by introducing the dust fields as dynamically coupled observers. It is only the 
total Hamiltonian and diffeomorphism constraints, which involve the dust as well as the gravitational and 
scalar field contributions, that is still vanishing. Moreover, as shown in [^, both H and Nj are conserved 
quantities because their Poisson brackets with the physical Hamiltonian equal zero. Now it is easy to see 
from Eqs. and (l5T]i that 


H = 


P(<I.) 


F\^)Qi,P^^ + 


Vdetg 


pjk _ ydet gp(<l>) ^^_iykmn 


2N 


N 


{Qr 


(4>-P-4>) 3(P'(4>))^-F2P(4>)P(4>) 


{C^Q)mn + p^^^ 


N 


(53) 

(54) 
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Taking the second order time derivative of Eqs. (1^1) and (I^T|) and substituting Eqs. (ISD]) . 
into them, we get the following equations: 


$ = 


N 


(^\/detQ^ 


^/dOTV N V 


\/det Q , 


N -^det Q \/det Q 
3(E'($))2 + 2F($)ii:($) N ’ 2 VdetQ3(F'(d>))2 + 2E($)K($) 


(i - + 2{M) + (£^3-) - {£^(£^3.)) 

C{t, a) 


N ' ' N ' 'TV 

3 F'($) 


+ 


+ 


2F($) 

3(F'(^>))2 + 2F(d>)K(d>) 


F'($) 

3(F'(^>))2 + 2F(d>)K(d>) 






(iVv/di^) 


- V det Q 


F(ch) 


(Or. - i^-Q)rs) [Qtu - {C^Q)tu 

-F'{^)Q^^{Q^n - {C^Q)mn) (i - 

+ { 2 N^K{^) + + 3iVV($) 

TSiV^D* AE(d>) + 3N{D^N)DiF{^) - N'^F{^)R^‘^'^ 
1 iV 


2 v/det Q 


:H{a)Q^^NmN„ 


Qjk 


N 


(\/det(5^ 


(55) 


.^ir ,vain7,N /. 

x/dOTV N ’) ^ 


N v/det Q \J det Q V N ') ' N 

2 K{^) 


FQ^^[Qm, - [Qnk " {C.-Q)n^ + 


+2JVD^D,V + 2(£ja) + (£ jQ),., - (£;;(£-<?)),., - (%<= - (£ 

2N 


N ' TV 

.rm,rn , 


Af 


TV ^ 


+ EW^ [iV(I^W)AF($) + iV^n‘AE(cl>) 

-2K(cl))F"(cl>) + i^^(cl))F-(c^) ■ ^ 2 K{<!>)Qj, 

3(F'($))2+ 2E(d>)K(d>) '' TV ^ VdiFQ3(F'($))2+ 2F($)it:($) ^ ’ 


N{D^N)DiF{^) + N^D^DiF{^) + A^V(^>) 

K{<^)Qjk 


2F'{^)Qjk 


’3(F'(^>))2+ 2F(^>)ii:(d>) V 




+K($)Q^''^ci>,J-^(iVv/d^)^ 


>(^>)[3(F'($))2 + 2F($)ii:($)] 


y/detQ 

^ (Or. - {C-Q)rs) {Qtu - {Cj^Q)tu 

-E'(<h)g-"(0„2n - iC^Q)mn) ($ " 

+ { 2 N'^K{<^) + + 3iVV($) 

+5iV2E>* AE(^>) + 3N{D^N)DiF{<^>) - N‘^F{^)R^^^ 
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( 56 ) 


+ -^^=H{a)Q^^NmNn , 

2VdiFQ ^ ' 

where C'(r, cr) is now expressed as a function depending on the configuration variables and their first order 
time derivatives given by 

C _ (C^Q)^s) {Qtu - {C^Q)tu) - v^d^F(cb)i?(3) + 

-^^^-F'mQ^^{Qrs - {C-Q)rs) (i - + 2v'ditQg^-^i),AF(ci>) 

+ \/detgy(«>). (57) 


So far we have derived the evolution equations for the gauge invariant 3-metric Qjk and the scalar field 
d>. Notice that if we choose in the equation of motions shown in (j55|] and (j56p F{^) = 1 and K{^) = 
all terms containing F'{^), F"{^), K'{^) naturally vanish and the remaining terms reduce to the evolution 
equations of the observables for a physical system consisting of dust, gravity and a minimally coupled scalar 
field. In this special case the evolution equations in (1551) and ([561) agree with the evolution equations for the 
observables shown in equations (4.23) and (4.24) in 


3 Linear Perturbations of STT on a General Background 

In this section we will use the gauge invariant variables {^,Qjk) and their equations of motions, obtained 
in last section, to derive the evolution equations for their corresponding linearly perturbed variables. For 
this purpose, we first split the configuration variables into a background and a perturbed part 

$ = Qjk = Qji^ + SQjk, (58) 

where and Qj^ denote the background variables satisfying the evolution equations shown in (f55]) and ([56]) 
and are thus solutions of the classical scalar-tensor equations. From now on we denote all occurring back¬ 
ground variables with a bar. Introducing the dust as dynamically coupled observers has the consequence, 
that the shift vector Nj and the lapse function N are no longer arbitrary but become fixed as particular 
functions on the reduced phase space, spanned by {Qjk, PF, <b, 11). Therefore N, Nj are not treated as inde¬ 
pendent variables and the same is true for the associated perturbations. Since we consider the perturbations 
of the second order evolution equations for Qjk and $, the lapse function N and the shift vector Nj can be 
understood as functions of Qjk,^ and their associated velocities. The explicit form of their perturbations 
reads 

_ n. _ 7] 

Nj = Nj +5Nj = -^ +5Nj, N = N + SN = = - ^6Qjk + =5Nj. (59) 

Note that in principle 5Nj can be expressed in terms of 6 Qjk and (5<I>, but we keep 5Nj here because as 
proved in both 5Nj and 5H are conserved quantities, while a conserved quantity is fixed at an initial 
time, it remains a constant during the evolution. Keeping 5Nj here helps to formulate the final evolutions 
equations for the linear perturbations in more compact form. Likewise we also keep the explicit expression 
5H in the following for the same reason. 

Since the evolution equations (|55]) and (l56]l are already very complicated, it is expected that the evolution 
equations for and 6 Qjk will be even more complicated. In order to express them in a concise form, we 

^Note that the definition of the potential V{(j>) used in differs by a factor of 2 and this needs to be considered in order 
to have an exact agreement between the equations. 
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notice that there are three important common terms in Eqs. (1551) and (1561) . namely C{T,a) and all terms 
inside the last two square brackets. We introduce the following abbreviations, which agree for both equations: 




N 


- V det Q 


(^N\J det , (60) 


Y = 


_ E'(<h) r^_li rstu 


{Qvs - {C-Q)rs) {Qtu - {C^Q)tu) - F'mQ'""{Qmn " {F^Q)mn) 


d> 

N 




+ {2N^K{<I>) + + 3XW($) + 5N^DWiF{^) + m{D^N)DiF{^) 


(61) 


The next step is to derive the equations of motions for the linear perturbations. For this purpose we use 
the following identities: 


[G-Y'" {Qts - {C^Q)rs) [Qtu - {C^Q)tu 


= 2 


(Qrs - {F-Q) 
V N 


G ^ 




■ 1 ] turj 

Gr 


Q'''(Qtu-iF-Q)tu 

V N 


+ [(? 


-2 


G 


rs]n 


(<3„ - (C-^Q) 

Q'""((3m„ - (£jQ)mn) (i - £-«) 


<?’X[ej_,+2e,„—(g'V 

Y‘[Qj„]j + 2Q,j^(q' 


^Qjk 


5Nm. 


Q''"Q’'"{Qmn - {F-Q)mn) 


-g”^"(¥-£^¥) 


N 




(g^„ - (£^g). 


^Qjk 


+ 


-g^^($-£-$) 

N 




'' dx^ 


-g'"[g™"^,n](g,fc-(^-Q).fc) 




+ 


g"”(g^„-(£^g)^.)(^-£ 


dr N' 


(5$, 


]^[DnF{^))DmQ"" 


6Qjk + [D,DkF'{<^)] 6^, 


6[DjDkFm = [-(D„F(cl>))22(fcg™"] <5g,)^ + 

NQF{DjN)DkF{^)\ = -ig™’"(;D„F(¥)):D^(]V']v'^) -ig^”^g"'^(;D^]v')(;D„F(¥)) 

T 


+ 


^Qjk 




+ 


(zi,w')g'"i?fcF'(ci>) 


5d>, 


N 




1 -/V ,yr;^n N N N-r77=7 “jTr^lTr^ 

(g H -ZZ2 — )HGjkrsF^ N + 


1 X 


2 ydetg 


N 


2 v^detg 


HN N Qjf^ 


(62) 


(63) 

(64) 

(65) 
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+2' 


N 


+ 


+ 


+ 


vdetQ 
1 N 


HGjkrsN N^Q 


5Q ran + 


- 2 - 


N 


TT^^t-rrSm 


vdetQ 


HN'N^Q^"^GrtsU ^Qk) 


HN"^N'^Q^^ 6Qjk 


2 y/ det Q 
N Jl 
\/det Q N' 

= GjkrsN N 


^HG.krsN N -2 


N 


y/detQ 


HGjkrsN^Q' 


6N„ 


y/ det Q 


6H, 


( 66 ) 


where the derivatives inside the square bracket act on all terms including the perturbed variables to its right 
side, while derivatives surrounded by the round bracket like (DiDj...) act only on the elements inside it. 
With these identities, the perturbation of the expression X in (|60l) can be expressed as 


6X = 




^jk 


d 1, 




N 


\/detQ 


(Wdet^. +K{^) 


dx^ 

^ (nJ^) Iq'"' ^ 


- -\/ det Q 


dx^ 


5<1> 


+ 


-Wn'^ 




+K(^) 


dx^ 
1 N d 


® (Q^”Q*“¥„) 


1 N —ik — ! -= 

+ ^^)[A'Vdetg],™ + — 

2 v^det Q N ^ 2 ^det Q dx 


N 




+ 


VdetQ 

2N^ (^K(^)lTDi'^ + i^(¥) [g 


^Qjk 


1 ,m 




N N'' 


- s/ det Q n‘^ 


[Ny det Q) 


6N„ 


(67) 


and the perturbation of the expression Y in (I6ip yields 
F'i^),T^-hrstuffs 


6Y = 


[G - {C^Q)rs) [Qtu - {G-Q)tu) " F(cl>)"(cl>)g™" - {C^Q^n) (' 




N 


-F' 


''(^)g"”(g„^n - {G^Q)mn){^ " Cy_) + iV'g™"$,^$,42iVV($) + iv'^if'($))' 
+2(2ii:(¥) + + 3]vV'(¥) + bx'^VDiF'i^) + m{D'N)D,F'(^) 


-N W ’F'{^) 




+ 


n^) (i^ 


{Qrs-{^^-Q)rs) [G ']~(^ -/:-)-[G Y^'"Q'\Qtu-C--Qtu) 
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-F'(c^) [ - - {C^QU) - c^<^) + - /:-‘i>)(|- - £-) 

V AT / N N OT N 

+0"" (4«; - {J^-Q)vw) [Q”'"¥ .iV"] 

-Q"-{^-C^^)[-Q-^iriQ^J,-Q^,^iQ-^ir^ 

-3v(^)7r"7r +5[ - ir'ir'^^iDyD^Fi^)) + 

-rT2-^mr ■r-i/“^\\“7T -T^nw -r;rrp‘ 2 '- 7 =^V-T:rnW-=r -rn/^w 

-N Q {D^F{<^))DrQ -NQ Q {D^Dy,F{<^)) 

+3|^--Q (D„F(3>))L>^(A^ N)--Q Q {DyN ){D^F{^)) 

+F(^)R^^^7rW - N^Fm\[G~^]^^^^DyDy, - IT''' 

Itt/ \T=7'^^'^A^^rT ITT r 1 ^ .-T^mn N N ,1 1 N —— 

+ 0 -^(^)g Nv^w ~ 7T— / _ (g H-=2—) ~ O— ; _ H{a)N N 5Qmn 

2 L 2 yMetg iV^ j 2 y^detQ 

+ -^{Qrs - ic-mrs)[G-Y^-[Q^\Qyj,t + 2g,„^(g”^*)) 

-f'(¥) - - /:^¥) \^\Qyj,r + g«,.^(g"’') + Q.r^(g™")] 

-g"" (4«, - [g”'"^,n] 

+2]v™(2i^(¥) + (¥))¥’*¥ i + 6y(¥)F™ 

+10[W"Q""(DyDy,Fm)] +3\^"(Dy,Fm)Dy'ir] - 2F(¥)72^"^]V" 

+ ^ ^^ Hia)'^"NrNs{^) + ^£=:ff(cT)g”^’^]vJ <5iV^ 

2 Vdetg N Vdetg 

+ [ ^ J 5//. (68) 

[2 v^d^ J 

Furthermore, we will need the perturbation of C{T,a) that is given by 

-JC = i [g"'] (g,, - {C^Q)rs) (Qtu - iC-Q)tu) - Vd^^^'^F'(¥) 

+yd^iF'(¥)g'%¥fc + 2yd^iF(¥)g'%^ - F'(¥)g”^"(g^, - {c^QUn) (^ - 
-F"(¥)g”^"(g,,, - (/:^g)^n) (^ - + ^^!§3iF'(¥)(¥ - £^¥)' 

V Af / N N N 

+^ V^tg ^(-) ^ 2Y^detgg^'^:D^-:DfcF'(¥) + 5^ 

+ [ig^'^G + [g-' 1(4, - (G^g)..) - {c^Q)tu) 

2 N \ AN L J V Af / V AT / 
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- {c^Q)rs) 


+ 


V det Q , 

—=2 —^ 

N 

\J det Q F{^) 


+ [G 


+ 


V^dst Q I^ \ I^ Y' ^^2 


TV iV 

9r TV 

n^nJ,^ + 2QinT^(Q'"^" 


N 


^sk . 


(q„ - (£-<3)„) [(G-'|’'*^'=(^ - £-) - [G-'|‘“''<3“(<5,„ - £-<5,J 


'dx'^- 


_ v^k(¥)¥''’?’‘ - f'(¥) [ - Q”^e''”(i5„„ - (CJ)U„) ( 


N 

Vk, 


- Y^det QF($) [[G Y^^'D^nDn - R’’"] 

t - /:-¥) 


N 


d 


N OT N 


+ Q"" {Qmn - {J^-Q)mn) [Q''’^,rN’^] 
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With the above formulas, we have all intermediate results available in order to derive the evolution equations 
for the linear perturbations (5<h and SQjk- We will start to discuss the case for 5<h. First we write the evolution 
equation for 5<1> in the following form 

+ B^^6Qjk + C^SNm + BdH. 

This is always possible with appropriate choices for the coefficients A, B-^^, C”^ and D. Using the explicit 
form of the coefficients A, B-^^, C™, D in the above equation we obtain 
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Notice that we did not write down the explicit expressions of 5C, SX and 5Y here in order to show the 
final equation of motion in more compact form. Now we consider the evolution equations for the linear 
perturbation SQjk- Analogue to the case of (5d> we can also formulate the final evolution equations in the 
following form; 


SQjk = U,fc,5^ + V6Qjk + WyySQmn + Xy^5Qk)m + YTk^Nm + Z5H, 

with appropriate choices for the coefficients \Jjk, V, W^”, X|^, and Z. Using the explicit form of these 
coefficients the final evolution equation is given by 
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and metric fields when the background metric are solutions of the STT equations of motion without any 
further symmetry requirements. Therefore these equations describe the evolution of the linear perturbations 
for general, and not only cosmological, spacetimes in the context of STT theories. The evolution equations 
in (j70p and (1711) are the generalization of the evolution equations (6.13) and (6.31) in to the framework 
of STT theories. However, for the purpose of practical use, in the next section we will choose a specific 
background, namely a cosmological one, and study the detailed behavior of the perturbation equations. 
This also allows to compare the results obtained in this manifestly gauge invariant framework with already 
existing results in the literature obtained in the standard framework of cosmological perturbation theory. 


4 Cosmological Perturbations on FRW Background 


In this section we specialize the general relativistic equations from the last section to homogenous and 
isotropic FRW backgrounds and compare them with the already existing results in the literature. For 
simplicity, we only focus on the spatially flat case. The equations are greatly simplified because the spatial 
derivatives of background variables in Eqs. (fTOll and m vanish, and we have Nj = ^ = 0 and N = = = 

1 + Q^^NjNk = 1. We denote Qjf. = A^5jk where A represents the gauge invariant scale factor. Then 
from Eqs. (I55p and (|56l) we obtain the following evolution equations for the background variables 
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e = H = A^ 


-6F($) 


- 6^F'($) $ + K{<^) ($)2 + E($) 

./I / ./I 


(74) 


Notice that e is a constant, because as mentioned above the physical Hamiltonian density H is conserved 
during the evolution. From Eqs. (I73p and (1741) we read the gauge invariant Eriedmann and Raychaudhuri 
equations as: 

(4)" = i(4+p). (75) 

2|-(4)2 = -i(P| + P). (76) 

where p® and P® denote the effective energy density and pressure of the scalar field, while p and P represent 
the correction terms caused by the non-vanishing physical Hamiltonian density. These quantities are defined 
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Given the evolution equations for the background variables we can consider their linear perturbations around 
an FRW background. The evolution equations for the linear perturbations have the following form 
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To study the detailed behavior of the perturbation equations (I77p and (I78D we will work in conformal time 
x^ = r] for which g^Q = —A"^. However due to the dust observers we always have grr = —N'^+Q^^NjNj, = —1 
where r denotes the physical time, meaning that we automatically work with the proper time of the dust. 
Using the relation dx^ = ^ we obtain 

9ti — — Qjk- ( 79 ) 

Following the procedure in standard cosmological perturbation theory (SCPT), we decompose the spacetime 
metric g^^, = g^,y+5g^i, into tensor, vector and scalar modes, parametrized by ten functions (j), B, E, if), Sj,Fj, hjk 
where the vector modes Sj and Fj are transversal with respect to the spatial Euclidian metric and the tensor 
mode hjk is transversal and traceless, 

9oo = (~1 + 2<^)H^, goj = [Bj + Sj)A‘^, gjk = 5jk + At^^i^djk + 2Ejfc + 2F(^jk) + hjk\- (80) 

Now from g^-r = —1, we immediately get ^ = 0 and therefore with the dust as dynamical observers we will 
always work in partly synchronous gaugqj- Furthermore since Nj can be expressed in terms of the other 
phase space variables, the functions B and Sj are not independent variables in our formalism. We remark 
that each mode in equations (I80p is already gauge invariant in our formalism, which is different from the 
case of SCPT. At the level of linear perturbation theory the individual modes decouple and we can study 
their evolution separately. Let us denote the derivative with respect to conformal time ^ by an upper prime 

'. Substituting the expressions (IHO]) into (1781) we obtain for the tensor modes 


, F’(<^)— / 

2'Hhjk + pi^'^ ^^hjk + hjk — Ahjk = 0, 


(81) 


where T-L A := 6^^DiDj. For the vector modes we consider the decomposition of the shift vector from 

(fHO]) given by 


6Nj = A[B,+Sj]. (82) 

Since 5Ni = 0 and in the case of the vector modes we also have H = 0, we obtain 

^Sj + A^ = 0 nSj + S'j=0. (83) 

dg dg j j 

^Here partly refers to the fact that in SCPT the gauge <p = B — Qis denoted by synchronous gauge, while longitudinal gauge 
means B = E = 0. 
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From (1781) we get for the vector mode the following equation 




2^(i,fc) + ^'FFij,k) + 2 
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(84) 


Denoting Vj = Sj — F- and using Eq. 


Eq. 
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Eor the scalar mode equation (f78l) leads to 
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where we used the field redefinition 
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( 86 ) 


(87) 


Eq. (1861) is of the form fjk + gdjk = 0 for appropriate definitions of the functions / and g. Applying the trace 
to this equation we have 3g + A/ = 0. Now operating with djdk on fjk + gdjk yields A{g + A/) = 0. Since 
due to our boundary conditions there are no zero modes of the Laplacian, we can conclude f = g = 0 and 
therefore the {.)jk part and 6jk part of Eq. ([86l) are independent of each other. Considering the separated 
contribution we get 
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( 88 ) 

djk 

(89) 


In SCPT the system does not involve any physical observers [^, the gauge invariant variables of STT are 
defined as: 

©A = -(5In -[n + (In ^ F(¥))'] [B - e') - {B - e')\ 

Qb = [^ + 5\n^/F(^] + [n + {\n^|F^))']{B-E), 

bZ = (5H+ h'(S-E'). (90) 

Notice that in Ref. [^, there is an additional term proportional to bN in the definition of 0yi. However in 
our model bE is not an independent variable and we have 


bN = - 
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Now for an FRW background we have = 1 and Nj = 0, and thus SN = 0, explaining why no contribution 
to 5N occurs in our equations. In the case of the scalar mode we have 5Nj = ABj and consequently 


dA 


B j + A 


dB. 


dr] ' dr] 


= 0 


T~LBj + B j — 0, 


(92) 


where the last relation follows from the fact that the Laplacian has no zero modes. Using Eq. (1921) . Eqs. 
and (l8^ are simplified to 

2[0a - Qsijk = 0, (93) 

Q'b +[n + (In /^)'] (0A + 29b)' + [2[n + (In /^)']' + [H + (In /^)']')0 a 


= -^[{E feA + E6Z' -F{E)A^{^^^)yz] 


(94) 


where )s d( Note that in a slightly different notation was used but the notations 




can be related to each other by specializing to E(<h) = 1 and identifying the variables and 'k in 27| with 
0A = 4> and 0^ = 'k. 

As the next step we derive the equation of motion for the perturbation (5<k. Using the field redefinition (1871) . 
the equation of motion for the scalar field perturbations in (I77p can be written in a more compact form as 




iVi. \/F(i))L - + 3{(ln V^) + l) (1“^’(2 ))h2 


V( 


A. 




(IS + 


- 4(1.1/F(=))k2-34 


1 1 1 ^ ■ 1 ,1^ 1 


C 


A-^^A' 2A2 

Using the definition (fMl) . it simplifies to 


2 F(H) v^detQ 


)• 


(95) 


dZ" + 2 [n + (In ^F{E))']6Z - A6Z + {^^^y^^F{E)AHZ + ^'0:4 - 2 {^^^y^F{E)A^e'A + 

^ ^ (96) 


= (-(^^)=4(^ -+ 4^^)=4(^-^7 - 42ri(i(^—)h ^ ). 

4V(“)^“A^ U 2^F{E)’-A^ ^ 4^F(H)^^VdetQ^ 


Note that compared to SCPT in the manifestly gauge invariant formalism we start with STT plus dust and 
thus increase the number of degrees of freedom by four. These four additional degrees of freedom are used 
to construct gauge invariant quantities associated with the spatial 3-metric qjk and the scalar field (f>, which 
we denoted by Qjk and ‘k respectively. In order to compare our framework with the results of SCPT we 
expressed 6Qjk in terms of the invariants 0 a,0b,Ij' and hjk, usually used in SCPT, and also introduced 
the field redefinition in Eq. (|87p to finally work with <5H and 6Z respectively. Now if we properly identify 
QA,QB,yj, hjk and 5Z with the analogue quantities in the SCPT formalism (see [27| for a detailed discus¬ 
sion), we realize that the evolution equations for the seven perturbed invariants 0 a, 0s, Vj,hjk and 6Z take 
a form almost identical to those equations that one obtains in the Lagrangian formalism using SCPT. The 
almost refers to the fact that in the manifestly gauge invariant formalism the dust as dynamically coupled 
observers has the effect that the physical Hamiltonian, which generates the evolution of the gauge invariant 
quantities 5Qjk and 5<k, is no longer a constraint as the case for SCPT. Instead it becomes a constant of mo¬ 
tion. This constant of motion reflects the influence of the dust and hence the dynamically coupled observers 
on the system. In case of Eq. (1961) these are exactly the terms on the right hand side of this equation, which 
are corrections caused by the non-vanishing Hamiltonian density e compared to the corresponding equation 
in SCPT. These corrections also encode the fact that in the manifestly gauge invariant model we have seven 
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physical degrees of freedom whereas for STT in SCPT the theory possesses only three physical degrees of 
freedom. This can be seen from the fact that in the case of SCPT the seven invariants 0^, ©s, Vj,hjk and bZ 
are still subject to constraints which reduce their independent number of physical degrees of freedom down 
to three. In contrast for the manifestly gauge invariant formalism the constraints have already been reduced 
by constructing the observables Qjk and <h. In particular in the limit of vanishing influence of the dust, 
meaning that its energy and momentum density vanishes, the conservation equations turn into constraint 
equations again and the corrections that occur in the evolution equations for the linear perturbations vanish. 
Then we obtain an exact match with the results obtained in SCPT. However the whole manifestly gauge 
invariant formalism is based on the idea that the influence of the dynamically coupled observers cannot 
totally be neglected and thus small corrections would always be needed to taken into account. As discussed 


m 


27| for the dust observer in the late universe these modifications decay as compared to the usual terms. 


As mentioned above the constraints in the SCPT become conservation equations in the manifestly gauge 
invariant framework due to the dynamically coupled observers since here only the total constraints consisting 
of STT contribution and the dust contribution vanish. First, we have two conserved quantities associated 
with the energy and momentum conservation respectively: 


ej{a) =-Cj{a), €{a) = H{a). 


(97) 


In the following we will derive the equations for the perturbed conserved charges bej and be around an FRW 
background dehned through bej = ej — Zj and be = e — Z where for FRW Zj = 0 and Z = H. We obtain for 
the linear perturbations 

bej = F{^)A{bQjk,k - bQkkj) - 2F(^)A{bQjk,k - bQkkj) - A^F'{mS^)j " F{^)A[AbNj - bNk,kj] 

- [2AA^F’(^) + A^F’A^W - 6A‘^Af’{^) + 2A^tK{^)] (5$)j, 

(98) 

and 


be = 


- 6F'($)A3(^)2 - 6F"($)A3^$ - 6A3^F'($)^ + A^K'{^){^f + 2A^K{<^)^-^ 




+2A^ F'mDDi + A'*—^ 

i^(^)J 


(5$ + 


^ z — B — ( AA 

[2A^F{A>) ^ ’ dr ^ ’ A 


+ 4F(¥)i(f™"—+ 2F'(¥)A¥(f”"^-^ bNm. 
ax"- ox^ 


^Qjk 


(99) 


As the next step we decompose the vector equation 

is the Green’s function of the Laplacian A, ai 
and (f82]) into ([M]) , the longitudinal part of Eq 


into the longitudinal part be''- = A ^bek,kj, where 
A“^ is the Green’s function of the Laplacian A, and the transversal part (fe-*- = bej — JeL Inserting Eqs. 

J 

gives 




( 100 ) 


and the transversal part gives 


APi = 


bei 


F($)A2 


( 101 ) 
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Written in terms of variables 0^, 0^ and 5Z, Eq. (|lU0p becomes 


QB+[n + {\n^jFm']QA + \^SZ 
Finally, Eq. (IMj) gives 


1 




F{<!>)A L A 




(102) 


A0b -3[n + (In V^)'] ©B + I [ - 


1 


4F(^>)A _ 


6e - e(<5 In v^E(4>) + 3(V’ + 6^F{<^)) + AE - 2{B -E')' + [n + (In y^F($))'] {B - e')) 


(103) 


Up to now we have derived the equations of motion for the background variables in (1721) . ([751) and ([76]) and 
the perturbed variables in (fHTl) . ([85]) . ([93]), ([Ml) . ([M|) and ([1011) . ([102p . ([103p in manifestly gauge invariant 
formalism in the case of STT. Aga,in we find if we choose E(<h) = 1 and K{^) = ^ the STT + dust reduces 
to GR + scalar field + dust in (27(], which allows us to test the correctness of our equations here. 


5 Gauge Invariant Hamiltonian Formulation in Einstein Frame 

It is well known that in SCPT of STT the Jordan frame is related to the Einstein frame by a conformal 
transformation and can therefore be transformed into each other by appropriate field redefinitions. However, 
it is not clear whether this kind of equivalence will still hold in the manifestly gauge invariant formalism 
since such a transformation need to be formulated at the level of the reduced phase space here. In the first 
part of this section we will prove this equivalence. In the second part, using this equivalence, we extend our 
results to slightly different reference systems. 


5.1 Conformal Equivalence 


Under the conformal transformation gab = E{(j))gab and the field redefinition ^ = yS] 
action (|42l) in Jordan frame is transformed into the following action in Einstein frame. 


E1M-]2 
FW ) 


+ 2 


K(4>) 

FW' 


the 


^ {^,g,dust) j ^ ^ V^l det (^) I 




— - d^x-\/ det(5) 


m i 


-g^'^UaU, + 


m 


(104) 


where gab is now taken as the basic variable, stands for the scalar curvature of the Einstein frame 
metric gab, 9°“^ is the inverse of gab, V is the covariant derivative compatible with gab, U(0 •= ^{4’iC)) and 
F{^) := E(</>(^)). Notice that the definition for the dust variables p and are those in Jordan frame. This 
means that the reference system remains unchanged. 

Hamiltonian analysis of action ([104[) is performed in appendix A, where we get the following physical 
Hamiltonian, 


H(r)= f d^aH{T,a), H(r, a) =a), 
J s 


(105) 
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with 


C, (r, a) = + I[DjE\ (r, a), 


C{T,a) = 


-jL={Q,rnQ,n - \Q^Qmn)P^^P^" - 
Met Q 


1 1 

+ 2 /-^ + 2 

A/detQ 


+ -VdetQ(g*^(Ad>)P,d>) + -^) {T,a). 


The elementary Poisson brackets read 


{P^Ha),Q^n{a')] = 5lji^5{a,a'), {n(a),H(M)} = 5{a,a'). 


The evolution of a general observable is given by 


= {H,0} 


= j^dV|^(c'(u){C(u),0}-Q^'"(u)C',(u){C,(u),0} + ^Q^-(a)Q"-(a)C,(a)Cfc(a^ 

= j^d^a{N{a){C{a),0} + N^{C,{alO} + {a){Q,,{a\0} + ^^{F(a), O}) (109) 

where the dynamical shift vector and lapse function are defined as 


PC, 

Ni = - 

^ H 


, FC ! - ~ - 

H ^ 


By calculating their Poisson brackets with H, we easily verify that Ij := —Cj{a) and e := H{a) still compose 
a pair of conserved quantities. Hence := also remains a constant during the evolution. 

Simply following the procedures in section two, we get the second order equations of motion for H and 
g^Tc as 


-jv (ydetg)'^ jv 
^ A/det g A/det( 


-)](=- C-S) + Q’%, [-^ (|JVv'det Ql„ - ; 


N N 


+iV^ [mB + - 5( (HV)] + 2(£j=) + (C.E) - J=)), 


N 


[pm 


N N 


-)) (4fc-(^- 


- {C-Q)mj) (pnk - (P^QU) + - 2Rjk) 

+ 2 ND,DkN + 2 {£^Q) +{£^Q) -{£^{£^Q)) 


NH{a) 

UetQF{E) 


~ ~ ~ r iV^ - -0 V(E) 1 

Gjt,„„V"*W" + (Jjt- ^=o+N \ J . 

^ 2 ,/dT 5 J 
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It is not difficult to check that Eqs. (11111) and (|112l) exactly reproduce Eqs. and ([5511 after substituting 
Qjk —>■ F{^)Qjk, —>■ •\/^( Thus the evolution equations in Einstein frame can be 

related to those in Jordan frame through the conformal transformation and field redefinition. Since the 
equivalence holds for general variables, obviously it also holds for the linear perturbed variables under the 
assumption that the transformations between the two frames are non-singular everywhere. Thus we conclude 
that the Jordan and Einstein frames are still equivalent to each other in the gauge invariant formulation. 


5.2 Generalizing to Different Reference Systems 

In above sections we chose the dust particles as the observers, whose equations of motion satisfy g^'^U^Uy = 

— 1. An interesting question is whether one can choose different observers, for example the ones satisfying 
g^'^U^Uy = —X{(j)) with an arbitrary (positive) function. To answer this question, we first generalize 

the action (llOJp to 

S {i-g,dust) = J det(5)| 

(113) 




1 


p 




r'^UyUy + 


m 


where J(^) and T(^) are arbitrary positive functions. Since it has been shown in the last subsection that 
the two frames are equivalent, we choose Einstein frame in which the equations look simpler. The physical 
Hamiltonian reads (see Appendix A) 

U{t)= f d^aH{T,a), R(r,a) = - Q^^CiCj{T,a), (114) 

J s 

where expressions for Cj and C are the same as Eqs. (jl06l) and (jl07l) . The dynamical shift vector and lapse 
function are defined by 


N, = - 

^ H 


N=^ = JL + Q^iNiNj. 
H ^ 


(115) 


In FRW background, we have Nj = 0 and N = L{’E). The equations of motion of background variables 

read 


. A ^ o 

3(^)^ 


3(- 


(In - 3(-)E - L{E){^^-^)L - 

1,-2 L{E)ViE) _ I 


-1^1^ L(H)1/(H) 3^2_ l(H)E_(E) 

4^2^ F2(H) 2^ F‘^{E) ^ 


AJL{E)A^ 


(116) 

(117) 

(118) 


where we denote Qjk = A'^djk and e = H. The reader can easily check that the above equations reproduce 


the results in Eqs. (|72|) . (1751) and (176)) by setting L{E) = F[E) and replacing A with yF(H)A. Eor the 
linear perturbation equations, using the following decomposition 

SJV = — / ^ dJVj = 'JL{X)A{Si + Bj), SQjk ■= A^[2'il)6jk + 2E jfc J- 2F(^j k) + (119) 


2JL{E) 
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the evolution equations of the tensor and vector modes in the conformal time frame ^ ^ ^ are 


‘I'hLhjp. -|- hjj^ — ^hj]^ — 0 , 

SJI 

i2 ’ 


= 2nVj + Vj^, = o, 


(120) 

( 121 ) 


where 71:=-^. The scalar mode contribution gives 


0^ + HQa + sz 




-- {B - E'j) 


( 122 ) 


AeB-me'B + ][- A 


V{E] 




1 


Se - l{6 In VF($) + 3^^ + - 2{B - e')' + n{B - e')) 


(123) 


4^LiE)A L 

0s + nOA + 20b) + {2n + n^)QA = [(H )20a + E6Z' - , (124) 


where 0^ := -6 In - n{B - E ) - {B - E) , Qb :=-ip+ n{B - E ), 6Z := 6E + E {B - E ). It is 

easy to see that the above definitions are the same as those in (|90ll when L{E) = T(<I>(H)). The evolution 
of the perturbed gravitational scalar field reads 


SZ” + 2nsz - A6Z + (I^y^^A^6Z + H 0:4 - + 3H 0s 

= (_ Aplki(^ _ _ Apki(B - E)' - "i£>t ^). 

2 L 2 (r.)A 2 L 2 (r.)A L(^) 2^L{r.) 


(125) 


Again we can easily match these equations with the ones in Jordan frame. We find that the different choice 
with g^'^U^Uu = leads to the change of dynamical lapse function N and the correction terms. 

Now we compare our manifestly gauge invariant cosmological perturbation theory (MGICPT) with the 
standard cosmological perturbation theory (SCPT). For the convenience of readers, we list all the equations 
in table 1. We denote the non gauge invariant background variables in SCPT with lowercase letters and 
the gauge invariant linearly perturbed variables with a hat. Since the Jordan and Einstein frames are 
equivalent in our manifestly gauge invariant Hamiltonian formalism as well as in the standard formalism, 

we write all equations in Einstein frame where they take a concise form. In Einstein frame A := F{E)A, 
a := a := E{^)a, h := definitions for 0^, 0_b and bZ are 

given in Eq. (1901) . From the comparison we see that these equations match each other precisely provided 
that the Hamiltonian density e goes to zero. 
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Table 1: Comparison between MGICPT and SCPT of STT in Einstein frame 



Manifestly Gauge Invariant Formalism 

Standard Formalism 

Background 

Variables 

s'- 2m' 

VF2(=)/= 2Fi(E)A 

QU +-m) ^2(g) 

Gh + 2(0 + « - 0, 

2h + m a? ^2^1^ = 0, 

Tensor-Mode 

+ ^jk — ^hjk = 0, 

2^^afe + ^ab — = 0, 

Vector-Mode 

Au = -if. 

+ Kk = 

ak = 0, 

2hVa + = 0, 

Scalar-Mode 

0J5 + HQa -h SZ 

L J J 

AOB-sne'^ + il-A^^eA 
m'sz'-A^(^y^sz] 

= ^ - (5e e( (Mn/Fp) 

AF^CE)A\_ V V / 

+3^^ + AE) + 2e{B - E')' - h{B - E') , 

0B + n{QA + 20b)' + {211' + n^)QA 
= -i [(=>0^ + E5Z' - AH^^y^5Z ], 

®B “f ^0A + — 0; 

A0B-3F0'B + i[-d2|^0A 

Aiz' 

= 0, 

®B h{&A + 2@b) + {2h + h?)QA 

= -i [(!>0a + !'<5Z' - a2(m)' , 

Perturbed 

Gravitational 

Scalar Field 

5Z" + ZHSZ ASZ + ( )L^A^SZ 

0^ 2(^2(=))s44 0^ + 3-0B 

= (i(54)k^(B - £'))' + i(^y-m(B - E')' 

Flsy'^(2(F(sy)s^/d^)- 

,5Z" + 2/i,5Z A,5Z + ( y^^d^6Z 

+? 0A 2( ^20 )^a ©A + ©B 

= 0. 
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6 Concluding Remarks 


In this paper we extended the manifestly gauge invariant Hamiltonian formalism introduced in [2^ to the 
case of STT. Our results are summarized as follows; First, we derived the Hamiltonian equations of motion 
in the manifestly gauge invariant framework using the Brown-Kuchf-dust as reference fields in the Jordan 
frame. Afterwards we used these equations to derive the evolution equations for the linear perturbations of 
the 3-metric and the gravitational scalar field on a general relativistic spacetime background. These are the 
generalization of the results obtained in to the case of STT. Then we applied our general result to the 
special case of a flat FRW background. This allowed to compare the results obtained in our formalism to 
the results from standard cosmological perturbation theory (SCPT). Likewise to the analysis in 2^ where 
gravity and a minimally coupled scalar field were considered, also in the case of the STT we have shown 
that our results and the one from SCPT exactly agree in the limit when the energy momentum and its 
perturbation of the dust vanishes. However, taking the idea of the manifestly gauge invariant framework 
seriously this limit is rather artificial since the dynamically coupled observer will always have an imprint on 
the system and here these are exactly the computed corrections compared to the results of SCPT. Second, 
we analyzed the question whether the Jordan frame and the Einstein frame are still equivalent to each 
other in the manifestly gauge invariant Hamiltonian formalism by showing their equations of motions can 
be related to each other through conformal transformations. Third, we generalized our results from the 
original reference system to a slightly different reference system in the Einstein frame and we showed that 
these equations are consistent with the ones in the Jordan frame. 


There are several possible directions to extend our results. First, since we have derived the linear per¬ 
turbation equations upon a general background, the direct extension is to choose a background different 
from the flat FRW case and study the cosmological perturbations in STT. For instance, the homogeneous 
but anisotropic Bianchi models are interesting candidates. Second, since STT can include a large variety 
of modified gravity models, we can choose some specific functions of F{^), K{^) and E(<h) to study the 
detailed physical predictions in the inflation period and compare them with the predictions of SCPT. Third, 
our formalism has the advantage over the procedure used in the context of standard cosmological pertur¬ 
bations theory that we can directly get the gauge invariant perturbed variables from the full metric to any 
order without worrying about the difficulty in constructing the higher order gauge invariant variables as in 
the standard way, thus it will be very convenient to use our formalism to study the non-linear effects such 
as the non-Gaussianity in the primordial perturbations. 
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A Lagrangian and Hamiltonian Analysis of the Action in Einstein 
Frame 


In this appendix, we will perform the Lagrangian and Hamiltonian analysis of STT in Einstein frame and 
derive its physical Hamiltonian. The action which we consider is the following, 


^ {^,g,dust) j ^ ^\/l I 




d^x-\/\ det( 5 )| 


J{i) V 




m 


(126) 
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1. Lagrangian Analysis of the Dust Action 

Recall that the dust velocity field is defined by and the energy-momentum tensor of 

the dust reads 


rpdust 2 5Sdust _ P 

~ vIditM ~ >^(0 


U U _L_ 

^ " 2J(0 




mr 


The variation of action (jl26p with respect to p gives 

r^UaUp + 


m 


= 0 , 


(127) 


(128) 


thus the energy-momentum tensor reduces to 


T — U TJ 


AO 


(129) 


which coincides with that of a perfect fluid with energy density and zero pressure. Similarly, the 
variations of action (I126p with respect to Wj, T and yields the following equations, 


CjiiiT = 




-L(0' 

With these three equations, the geodesic equation in Einstein frame is modified to 


U = 


(130) 


(131) 


thus unlike in Jordan frame, in Einstein frame T no longer defines the proper time along the dust flow lines 
and there is a correction term proportional to the gradient of ^ in the geodesic equation. The conservation 
of the energy momentum tensor transforms as 


The total energy momentum tensor of the scalar field and dust reads 
2 


(132) 


rptot . _ _ 

^ flU * 




The variations of action (jl26p with respect to the scalar ^ gives 


2J(0h2(^)^’^ 


Substituting Eq. (11321) into 


we get 




WT« = - (D|), - = 0, 


(134) 


(135) 


(136) 


Hence the total energy momentum tensor is conserved in Einstein frame. 
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2. Hamiltonian Analysis 

Although the Hamiltonian analysis for the action (jl26p with J(^) = L{^) = 1 was already done in Ref. 
[iH ]. we draw a word sketch of the Hamiltonian analysis for the generalized case when J(^) and L{^) are 
arbitrary functions. All the notations remain the same as those in [j^ . 

By 3+1 decomposition, the dust action can be written as 

Sdust = J J dx^Vdet{q)njl^{- + q’^’^UaUb + J^)- ( 137 ) 

The conjugate momentum of T and read 

P- = ^ = (138) 

^ (139) 

equations ()138l) and (I139p impose the primary constraint 

Zj := Pj + WjP = 0. (140) 


There are also four other primary constraints, 


. 7 SSdust p, „• yj dSdust „ n n 

I := Z := — - = 0, P := := =0, p ■= z := — = 0, Pa := ■= ^-^ = 0. 141) 

op,t oWj^t ont dre“ 




The resulted Hamiltonian constraint of the coupled system of action (jl26ll is given by 


-tot xQeo I -scalar , -dust 

c = -\- c + c , 


where 


^eo ^ 


-scalar _ 


-dust _ 


,,, [qacqbd - ]:qabqcd]p''’'p‘"^ - \/det{q)R^^\ 
Vdet(g) 4 

^ [7s® + 


(142) 

(143) 

(144) 

(145) 


with Ua = -T^a + WjS^^a- Note that {p'^'^Aab, }, and {P,T} are three pairs of conjugate variables. 

The diffeomorphism constraint is given by 


-tot _ _|_ -scalar i -dust 

^a ^a ' ^a ' ’ 

where 

cr^-2qacDbP^^, = =-PH,. 

Now the total Hamiltonian is defined as 


H = 


J dx^ {p^Zj + pZ + pjZ^ + lAz + p°'Za + nc^°*' + rPc^a ^), 


(146) 


(147) 


(148) 
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where y, /i“ are Lagrangian multipliers. Now we consider the condition that the primary constraints 

should be maintained during evolution. Explicitly we have 


-tot 


z,t = {H,p} = -c' 


Za,t 


n 




Vdet(^ 


p^^Aei{q) J{i) 






=: c, 


Z,-t = {H,Z^}=pjP-n^PWj,a + n^^it^UbWj,a. 


(149) 

(150) 

(151) 

(152) 

(153) 


We can solve the constraints (|152p and (11531) by fixing the Lagrangian multipliers and pj. The new 
constraint (I15ip is of second-class, by asking the conservation of this constraint, we can fix the Lagrangian 
multiplier p. It is also easy to check that both and are conserved during the evolution. After solving 
all the second-class constraints, we get 


Wj = 

J p 


/i = 0, 1 = 0, p=^^M=i 

\/det(g) 






1 


m 


Finally, in the reduced phase space with 2 ; = 2 ^ = 0, we are left with only two constraints 

-dust 


-tot t^Q&o I -scalar , -dust 

c = c + c , 


-tot _ 1 -scalar 

I 


(154) 


(155) 


where 


-dust 


= pJrWb + 


1 


-dust ^ PT, + PjSi = -PUa. 

Substituting Eq. (jl57p into (jl56p and using 7^“®* = -|- c®“^“^), we can rewrite the constraints as 

-tot 


(156) 

(157) 


= P + h, 


h = y/L{^)\/^ - q°-^CaCb, 


cf = P,+hj, hj = -hS^Ta + S^Ca, 


(158) 

(159) 


where c := ■.= ^a° + and = 6^. 
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